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ABSTRACT 

Shock breakout is the brightest radiative phenomenon in a supernova (SN) but is difficult to be 
observed owing to the short duration and X-ray /ultraviolet (UV)-peaked spectra. After the first 
observation from the rising phase reported in 2008, its observability at high redshift is attracting 
enormous attention. We perform multigroup radiation hydrodynamics calculations of explosions for 
evolutionary presupernova models with various main-sequence masses Mmsi metallicities Z, and ex- 
plosion energies E. We present multicolor light curves of shock breakout in Type II plateau SNe, being 
the most frequent core-collapse SNe, and predict apparent multicolor light curves of shock breakout 
at various redshifts z. We derive the observable SN rate and reachable redshift as functions of filter x 
and limiting magnitude mx^um by taking into account an initial mass function, cosmic star formation 
history, intergalactic absorption, and host galaxy extinction. We propose a realistic survey strategy 
optimized for shock breakout. For example, the f;'-band observable SN rate for TOg'jim = 27.5 mag is 
3.3 SNe degree"^ day~^ and a half of them locates at z > 1.2. It is clear that the shock breakout is 
a beneficial clue to probe high-z core-collapse SNe. We also establish ways to identify shock breakout 
and constrain SN properties from the observations of shock breakout, brightness, time scale, and color. 
We emphasize that the multicolor observations in blue optical bands with ~ hour intervals, preferably 
over > 2 continuous nights, are essential to efficiently detect, identify, and interpret shock breakout. 
Subject headings: shock waves — radiative transfer — supernovae: general — stars: evolution — 
surveys 



1. INTRODUCTION 

In contrast to Type la supernovae (SNe), systematic 
observational studies of core-collapse supernovae (CC- 
SNe) have been restricted at a redshift z < 1 (e.g., 
iPoznanski et al.ll2007l ). This is because CCSNc a re typ- 
ically fainter than Type la SNe (e.g., iRichardson et al.l 
120021 ) ■ except for rare energetic supern ovae (hypernovae, 
e.g., SN 1998bw. iGalama et al.lll998D or recent lv-found 
extremely bright SNe fe.g.. lOuimbv et alll2009[ ). 

Recent improvements of telescopes/instruments and 
well-organized survey strategies make it possible to de- 
tect unusual core-collapse events at high redshift: for 
example, bright T ype Iln SNe (SNe Iln) at z > 2 
(jCooke et al. 120091) and ganima-ray bursts (GRBs) u p 
to z ~ 8.2 (jSalvaterra et al.ll2009HTanvir et alll2009[ ). 
With the use of the high-z accessibility of these events, 
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a star formation history (SFH) and an initial mass 
funct ion (IMF) at high redshift are intensively studied 
(e.g.JKistler et al.l [2009: Wang & Dai 2009; Cooke et aU 
|2009( ). However, special conditions are required to re- 
alize such event s, i.e., a dense circu mstellar matter for 
an SN Iln (e.g.. iChugai et al"1l2004D and a fast- rotating 
progenitor for a GRB fe.g.. lWooslevlfT993l) . Thus they 
can not be the main constituents of CCSNe and the SFH 
and IMF estimated with them could involve large biases. 
Hence, ways to directly detect normal CCSNe at high 
redshift are required to investigate the nature of major- 
ity of CCSNe and estimate the SFH and IMF with smaU 
biases. 

The bolometrically-brightest phenomenon in the SN 
with a Shockwave is shock breakout. In a CCSN explo- 
sion, an outward Shockwave forms around a central rem- 
nant by depositing released gravitational energy. The 
Shockwave propagates through a stellar envelope to heat 
it up and accelerates its expansion. Since the star is 
optically thick, the Shockwave cannot be electromagnet- 
ically observed until its emergence from a stellar surface. 
When the Shockwave approaches the stellar surface at a 
distance with an optical depth < 10, radiation from the 
shock front starts to leak out and a hot fire ball sud- 
denly appears to emit a bright soft X-ray and ultravi- 
olet (UV) flash with a quasi-blackbody spectrum {T > 
10^ K). This phenomenon is shock b reakout which has 
been theoretically predicted by, e.g., iKlein fc Chevalieil 
(|1978l ). Its duration strongly depends on the presuper- 
nova radius. Brightness rises in several seconds to sev- 
eral hours and declines in several ten seconds to sev- 
eral days. Theoretical studies have suggested that the 
peak bolometric luminosity exceeds 10*^ erg (e.g., 
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iBlinnikov et al.ir200Cl[ ) and that shock break out is observ- 
able even if it takes place at z > 1 (e.g., iChugai et al.l 
120001) . However, the short duration and soft X-ray/UV- 
peaked spectra had made detection of shock breakout 
difficuh for a long while, except for fortunate detection 
of its tail in nearby SNe in U band (e.g., SN 1987A 
ICatchpole efanirOSTl: SN 1993J. IRichmond et al.l ll99l 



SN 1999ex. lStritzinger etalW^K ' 

First detection of shock breakout in the rising phase 
was obtained serendipitously in 2004-2008 and reported 
in 2008: Type lb SN 2008 D in NGC 2770 (di stance 
d = 27 Mpc, [S odcrbcrg et al.l 120 08: Mazzal i et al.l 
[2008t iModiaz et al.li2009; .Malesani et al.ii2009 '). Type II 
plateau SN (SN II-P) SNLS-04D2dc (reds hift z = 0.185, 
iSchawinskijFaDIlO^LGczari et al."2008'), and SN II-P 
SNLS-06Dljd {z = 0.324, Gczari et al. 200^). 

For SN 2008D, an X-ray LC from br eakout and opti- 
cal LCs after a tail were o bserved fe.g.. ISoderberg et al.l 
120081: IModiaz et all [20091 ). The X-ray LC rises in - 
60 sec and declines in ~ 130 sec but it is not clear whether 
the X-ray spectral energy distribution (SED) is thermal 
or nonthermal. For the SNLS SNe II-P, UV LCs of shock 
breakout and optical LCs of plateau were observed (e.g., 
ISchawinski e'Fani200l[G"ezari et al.l[2008l) . The UV flash 
of SNLS-04D2dc rises and declines in several hours and 
the subsequent UV LC shows rebrightening in several 
days. However, the UV data have too low signal-to-noise 
ratios to obtain SED.^ 

The observational papers and subsequent papers 
present theoretical models for shock breakou t; an- 
alytic models for SN 2008D (e .g., [Spderbe rg et al.l 

[200l IChevalier fc Fran sson 200l IModiaz et al.l 12001) i° 
and hydrodynamics calc ulatio n with two-t emperature 
radiative diffusion (jSchawinski et al.l I2008D and one- 
temperature radiation hydrodynamics calculation 
coupled with non-local thermodynamic equilibrium 
(non-LTE) spectra l c alculation for S NLS-0 4D2dc 
(ICezari et al.l I2008D . iTominaga etaLl (|2009t ) have 
performed a multigroup radiation hydrodynamical 
calculation for the first time and successfully con- 
structed a self-consistent radiation hydrodynamical 
model for SNLS-04D2dc. They have presented that an 
SN explosion of a 20AfQ star with an explosion energy 
1.2 X 10^^ erg reproduces well the UV-optical LCs of 
shock breakout and plateau. Also, the successful model 
demonstrates that an SN similar to SNLS-04D2dc is 
detectable at z = 1 with 8m-class optical telescopes. 

The UV-bright shock breakout is followed by a plateau 
phase. It appears as an SN II-P that is most frequent 
among core-collapse SNe (e.g., Mannucci et al. 2008; 
lArcavi et alllMot llJetaII[20ia ISmith eTalKSOlO ). In 
contrast to an SN Iln and a GRB, the formation of a 
plateau phase does not require any specified conditions 
other than the presence of a thick H envelope. There- 
fore, the SFH estimated with shock breakout in SNe II-P 
should have smaller uncertainties and biases than those 
with SNe Iln or GRBs. Moreover, observational prop- 
erties of shock breakout, the brightness, color, and time 

^ The UV observations of SNLS-06Dljd have too poor signal-to- 
noise ratios to extract the characteristics of shock breakout. 

The subsequent LCs and spectral evolution of SN 2008D are 
reproduced by the radiative tra nsfer calculation wit h homologous 
hydrodynamics evolution (e.g., I'l'anaka et al.|[2009ai rEIV 
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Fig. 1. — Presupernova density structures of (a) Z = 0.02 stars 
with Mms = 13M0 (red), Mms = 15Mq (green), Mms = ISA/© 
(blue), Mms = 20Mq {magenta), Mms = 25Mq (cyan), Mms = 
30Mq (black), and Mms = 4OM0 (gray), and (b) 2OM0 models 
with Z = 0.001 [red), Z = 0.004 {green), Z = 0.02 {blue), and 
Z = 0.05 {magenta). 



TABLE 1 

PaOGENITOR MODELS. 



Mms 


Z 


A/prcSN 


-RproSN 


-J^preSN 


[Mq\ 




[Mq] 




[lO^Lo] 


13 


0.02 


12.7 


564 


5.57 


15 


0.02 


14.1 


507 


4.69 


18 


0.02 


16.7 


713 


8.81 


20 


0.02 


18.4 


795 


11.1 


25 


0.02 


21.7 


1200 


23.9 


30 


0.02 


25.0 


1360 


31.4 


40 


0.02 


21.7 


1660 


53.8 


20 


0.001 


19.7 


756 


18.6 


20 


0.004 


19.5 


751 


13.4 


20 


0.05 


17.3 


1050 


15.7 



scale, depend on properties of the SN and its progenitor, 
the explosion energy E, presupernova radius i? nreSNi and 
ejecta mass Moj (e.g..lMatzner fc McKe3ll999[ ). Thus, it 
is possible to derive detailed properties of SN explosions 
from the observations of shock breakout and constrain 
an IMF precisely. 

Shock breakout is coming under the spotlight to probe 
high-z CCSNe but the observable properties are poorly 
understood. Therefore, in order to execute a shock 
breakout survey effectively, it is required to provide theo- 
retical predictions for observable quantities and propose 
strategies of observations and analysis to detect, identify, 
and interpret them. Hence, we perform multigroup ra- 
diation hydrodynamics calculations of shock breakout in 
SNe II-P with various Mmsj Z, and E w ith a multigroup 
radiation h ydrod ynamics code Stella (jBlinnikov et al.l 
1998, 2000, 12006D and present theoretical predictions of 
apparent multicolor LCs at various redshifts. Based on 
the theoretical models, we estimate the number of detec- 
tion and reachable redshift, clarify requirements on sur- 
vey strategies, and develop ways to identify shock break- 
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Fig. 2. — Spectral energy distributions at (a) t = 0, (b) t = 0.1 days, (c) t = 0.3 days, (d) 4 = 1 days, (e) t = 10 days, and (f) t = 100 days 
for different opacity prescriptions {red: original, green: including excited levels in bound-free absorption, and blue including inner-shell 
photo-ionization) . 

out and to derive the SN properties from the observa- 
tional quantities of shock breakout. 

In § [21 the apphed models and methods are briefly 
described. In § El results are sho wn. We present the 
multicolor LCs of shock breakout (S IB.ip and predictions 
of apparent multicolor LCs of shock breakout (§ 13. 2|) . In 
^ I3.3l we offer future prospects on shock breakout surveys: 
an expected number of detection and reachable redshift 

13.3. ip . dependencies on extinction and SFH I3.3.2p . 
requirements on survey strategies (§ l3.3.3p . ways to iden- 
tify shock breakout (§ 13.3. 4|) . and ways to constrain SN 
properties (i) l3.3.5p . In §|4l the conclusion and discussion 
are presented. 

2. MODELS & METHODS 

2.1. Progenitor model & explosive nucleosynthesis 

We adopt progenitor models with various Mus 
(= 13, 15, 18, 20, 25, 30, and 40Mq) and Z 
(= 0.001, 0.004 and 0.02) which are taken from 
lUmeda fc Nomotol (|2005| ). The stellar evolution calcu- 
lations include a mass loss depending on met allicity Z 
which is assumed to be proportional to Z^-^ fK udritzkil 
[2OOOI) . Since shock breakout arises at a thin surface 
layer with an optical depth r < 10, we adopt the stel- 
lar surface sufhciently outside which is as shallow as 
r = 0.001.^^ The density structures of the progenitor 
models are shown in Figures [T^ and [T|d. The proper- 
ties of progenitor models, A/mSi presupernova mass 
-MproSN, presupernova radius i?prcSN, and presupernova 

In order to justify the shallowness of the outer boundary, 
we perform a radiation hydrodynamic calculation for an envelope 
model down to t = 10~^ and confirm the consistency with the 
envelope model down to r = 0.001. On the other hand, we find 
that an envelope model down to t = 0.01 shows a bluer SED than 
that down to r = 0.001. This indicates that t = 0.01 is not shallow 
enough to treat SEDs of shock breakout. 



luminosity iprcSN, are summarized in Tablejl]-'^^ -RprcSN, 
being the most important ingredient for shock breakout, 
increases roughly monotonically with Z and MmSj ex- 
cept for the model with Mms — 15AfQ. The progenitor 
models are red supergiants with i?proSN = 500 — ITOOi?© 
producing SNe II-P. We note that the density inversion 
at the outermost layer corresponds to a super-adiabatic 
layer, where the temperature gradient is steeper than the 
adiabatic case. 

We calculate explosive nucleosynthesis adopting the 
same method as described in Tominaga et al. (200'^; the 
explosion is initiated as a thermal bomb, hydrodynam- 
ics is calculated including nuclear energy generation with 
the a-network, and a nucleosynthesis calculation is per- 
formed as a postprocessing. Since explosive nucleosyn- 
thesis ceases well before the shock emergence from the 
stellar surface, we perform a radiation hydrodynamics 
calculation of an explosion for a model with the abun- 
dance distribution after explosive nucleosynthesis and 
the hydrodynamical structure of the progenitor model. 

2.2. Radiation hydrodynamics 

We use the multigroup radiation hy drody namics code 
STELLA (jBlinnikov et al.l [199 8, 20o4 I2006D . The de- 
tail of STELLA is d e scribed in archival literatures 
(jBlinnikov et al.l [l998l I2OOOI [2001) and its reliability 
has been carefully verified by comparisons with analytic 
solut ions ([Matzner fc McKeel[l99 9^ 'Rabi nak fc WaxmanI 
2011), other numerical codes (Blinnik ov et al. I199S 



2003), and multico lor SN ob servations ( Blinnikov et al.l 
,1998, 2000, 2006; Chugai et a l. 2004). Here, we briefly 

We note that the progenitor structures, i.e., the relations 
among Mms, Mp^eSN, -RpreSN, and LpreSNi depend on the treat- 
ment of physics, e.g., ro tation, pulsation, mass loss, mixing length, 
and o vershooting (e.g., iLimongi fc Chieffill2006l : IVoon fc Cantiellol 
[20TOh . 
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Fig. 3. — Absolute bolometric (red) and multicolor LCs {green: FUV, blue: NUV, magenta: u' , cyan: g' 
the models with (a) Mms = 13Mq, Z = 0.02, E51 = 1, (b) Mms = 15Mq, Z = 0.02, E51 = 1, (c) Mms = 
(d) Mms = 2OM0, Z = 0.02, E51 = 1, (e) Mms = 25Mq, Z = 0.02, E51 = 1, (f) Mms = SOM©, Z = 0.02, 



Z = 0.02, E51 = 1 
(k) Mms = 2OM0 



(h) Mms = 
Z = 0.001, 



: 25Mq, Z = 0.02, Ssi = 4. 
E51 = 1, (1) Mms = 2OM0 



(i) Mms = 
Z = 0.004, 



25Mq . 
E51 = 



Z = 0.02, B51 
1, and (m) Mms 



10, (j) Mms 
2OM0, Z = 



yellow: r' , and black: J) of 
I8M0, Z = 0.02, £51 = 1, 
£51 = 1, (g) Mms = 4OM0, 



25M0 , Z ■- 
0.05, £51 



0.02, £51 = 20, 
1. 



describe the assumptions and procedures applied in 
STELLA and the setup adopted in this paper. 

STELLA solves the time-dependent equations implicitly 
for the angular moments of intensity averaged over fixed 
frequency bands and computes variable Eddington fac- 
tors that fully take into account scattering and redshifts 
for each frequency group in each mass zone. The 7-ray 
transfer is calculated using a one-group approximation 
for the nonlocal deposition of the energy of radioactive 
nuclei; we follow Swartz et al. (1995; see also iJeffervl 
[19981) and use a purely absorptive opacity for 7-ray. It is 
worthy to note that the 7-ray transfer does not influence 
the results in this paper because of no contribution to 
shock breakout from the radioactive decays. In the equa- 
tion of state, LTE ionizations and recombinations are 
taken into account. The effect of line opacity is treated 
as an expansion opacity accordin g to the prescription of 
Eastman & Pinto (1993; see also iBlinnikov et al.|[T998l ). 



We adopt 100 frequency bins dividing logarithmically 
from = 6 X 10" Hz (A = 5 x 10"^ A) to 3 x lO^^ Hz (1 A); 
such a number of frequency bins are enough to solve non- 
equilibrium continuum radiation and treat any SEDs ac- 
curately. We emphasize that there is no need to ascribe 
any temperature to the radiation. The coupling of multi- 
group radiation transfer with hydrodynamics enables us 
to obtain the color temperature in a self-consistent cal- 
culation, i.e., a luminosity- weighted blackbody fitting of 
SED. 

All previous computations with Stella employed 
the assumptions used in the code eddington 
([Eastman fc Pintd Il993f ) for bound-free transitions, 
in which all atoms and ions, except for hydrogen, are 
in ground states. Since new opacity tables for Stella 
will be released including excited levels in bound-free 
absorption (E. Sorokina in prep.) and inner-shell photo- 
ionization (P. Baklanov in prep.), we briefly examine 
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Fig. 4.— Spectral energy distributions at t = of the models with (a) Mms = 13Mq, Z = 0.02, E51 = 1 (red), Mms = 15Mq, Z = 0.02, 
E51 = 1 (green), Mms = ISM©, Z = 0.02, E51 = 1 (blue), Mms = 20Mq, Z = 0.02, E51 = 1 (magenta), Mms = 25M0, Z = 0.02, 
£51 = 1 (cyan), Mms = SOM©, Z = 0.02, E51 = 1 (black), and Mms = 4OM0, Z = 0.02, E51 = 1 (gray), (b) Mms = 25M0, Z = 0.02, 
E51 = 1 (red), Mms = 25M0, Z = 0.02, E51 = 4 (^reen), Mms = 25M0, Z = 0.02, E51 = 10 (f»/«e), and Mms = 25M0, Z = 0.02, 
E51 = 20 (magenta), (c) Mms = 2OM0, Z = 0.001, £51 = 1 (red), Mms = 2OM0, Z = 0.004, £51 = 1 (green), Mms = 2OM0, Z = 0.02, 
iSji = 1 (blue), and Mms = 2OM0, Z = 0.05, Esi = 1 (magenta). 
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5. — Color temperature evolution of the models. The colors and panels are the same as in Figure [4] 
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these effects on SEDs for a model with Mms = SOM©, 
Z = 0.02, and explosion energy E^i = i?/(10^^ erg) = 1. 
While the inner-shell photo-i onization cross-sect i ons are 
base d on fornmlae derived bylVerner et al.l ()1993l Il996al) 
and IVerner &: YakovlevI (|1995D ^'^ as in old eddington 
and STELLA routines, the case with excited levels treats 
all bound-free transitions, also for ground-levels, by dif- 
ferent fitting formulae as in code wmbasic^'* (Pauldrach 
[Ml see also E. Sorokina in prep.). 

The SEDs with three different opacity prescriptions are 
shown in Figure [5] Although there are small differences 
as < 0.3 mag in several frequency bins, the differences 
are diluted and diminished when SEDs are convolved 
with broad-band filters. This illustrates that Stella re- 
sults are robust with respect to different approximations 
for bound-free transitions at least for shock breakout in 
an SN II-P in which temperature is not extremely high 
{T < 10^ K ). Therefore, in th i s pape r we adopt a pro- 
cedure from lEastman fc Pintol (|1993f) . The opacity ta- 
ble in clude s 1.5 x 10^ spectral li nes from lKurucz fc Belli 
(fT995h and IVerner eFall (|1996b[) . 

3. RESULTS 

3.1. Synthetic multigroup light curves 

We perform the multigroup radiation hydrodynam- 
ics calculations for models with various Mms, E, and 
Z. Here, for simplicity, we set Moj to yield a canoni- 

See http://www.pa.uky.edu/~veriier/photo.html for details. 



http : //www . usm . uni-muenchen . de /people/ adi/Programs/Programs '. html" 



cal amount of ^^Ni without mixing to the envelope [the 
ejected ^^Ni mass M{^^m) = O.OTM©]. We note that 
the assumption does not affect the results because the 
radioactive decays do not contribute to shock breakout. 
In this paper, we adopt the AB magnitude system and 
the following filters: far UV (FUV) and near UV (NUV ) 
bands for G ALEX satcWitc (Morrissev ct al. '2005', 2007^, 
u' , g' , r', i' and z' bands (Fukugita et al. 1996), J, 
H, and K bands ([Tokunaga et al.n2002D . and F322W 
band for The Jam es Webb Space Telescope {JWST, 
IGardner et alll20(ffil) . 

The multigroup radiation hydrodynamics calculation 
provides wavelength- and time-dependent fluxes at an 
SN surface. Since lights radiated from different parts of 
the SN surface arrive at an observe r in a given direc- 
tion at different time (for details, see JKI ein fc Chevalieil 
p78"; 'Ims hennik eraP IT98ll : 1 nsman fc Burrowsl 1199? 
Blinnikov et al. 1 120021 120031) . a light travel time correc- 
tion^^ and a limb darkening correction in the Eddington 
approximation are applied. 

Absolute LCs of the corrected models are shown in Fig- 
ures [3^-[3hi. Figures [3^-13^, [3h-|3|, and [3)j-|3}ii show the 
LCs of the models with different Mms, different E, and 
different Z, respectively. The models with larger Mms 
and higher Z have larger i?prcSN, except for the model 
with Mms = 15Mq, thus having broader and slightly 

The light tra vel time effect in an asph erical explosion is in- 
vestigated in, e.g., ICouch et al] II2009I ) and ISuzuki &z Shigevarnal 
(201 Oil. 
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TABLE 2 

Properties of shock breakout. 



A/ms 
[Mq] 


Z 


E 

[10^1 erg] 


[McA 


'^c,pcak 

[105 K] 


^Imag -^rad.lmag 

[10-2 days] [10**erg] 
without corrections 


^pcak 

[lO^^erg s-l] 


'^c,pcak 

[10^ K] 


^Imag -^rad.lmag 

[10-2 days] [lO^Serg] 
with corrections 


^pcak 

[lO'l^'crg s-i] 


13 


0.02 


1 


11.2 


3.56 


0.651 


0.691 


10.3 


3.35 


1.09 


0.815 


6.99 


15 


0.02 


1 


12.7 


3.68 


0.512 


0.536 


10.2 


3.48 


0.915 


0.643 


6.64 


18 


0.02 


1 


15.2 


3.03 


0.994 


0.948 


9.41 


2.87 


1.46 


1.09 


6.90 


20 


0.02 


1 


16.8 


2.80 


1.19 


1.11 


9.11 


2.70 


1.68 


1.26 


6.91 


25 


0.02 


1 


19.9 


2.17 


2.94 


2.07 


6.78 


2.04 


3.40 


2.23 


5.99 


30 


0.02 


1 


23.0 


1.99 


3.51 


2.39 


6.52 


1.87 


3.97 


2.54 


5.86 


40 


0.02 


1 


19.6 


2.01 


5.16 


4.57 


8.55 


1.92 


5.60 


4.82 


7.90 


25 


0.02 


4 


19.8 


3.23 


1.27 


4.20 


32.5 


3.08 


2.28 


5.05 


21.1 


25 


0.02 


10 
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Fig. 6. — Comparisons among the models with corrections {filled 
circles), the models without corrections {open circles), and semi- 
analytic solutions (original: cyan shaded region and reduced: gray 
shaded region). The color of symbols represents Z = 0.02 models 
{red), Mms = 2OM0 models with Z = 0.001, 0.004, and 0.05 
{blue), and Mms = 25M0 models with E51 = 1, 4, 10, and 20 
{green). 

fainter peak. On the other hand, the models with higher 
E have brighter and narrower peak. While the bolomet- 
ric luminosities decline monotonically with time due to 



an adiabatic cooling, there could be rebrightening in ho- 
mochromatic LCs due to the shift of peak wavelength 
with time. 

The SEDs at the bolometric peak [t — Q, hereafter 
we set t = at the bolometric peak of each model) are 
shown in Figures H^-IH:. The SEDs ai t — Q peak in 
UV (A ^ 40 — lOOA) and have similar spectral slopes 
at A > 400A, while the luminosities at A > 400A are 
higher for larger i?prcSN and slightly higher for higher E. 
Evolution of color temperature is shown in Figures [5^- 
[S};. The color temperatures range from Tc ^ 2 x 10^ to 
Tc 5 X 10*^ K at t = depending on Mms and E. The 
SEDs and color temperature evolution depend on Mms 
(i.e., i?proSN) and E, while their dependencies on Z are 
small. 

The sem i-analytic solutions for shock breakout by 
iMatzner fc JVIcKea (|1999, ) provide radiation temperature 
?MM99, outburst energy i?MM99, timescale tMM99, and 
luminosity Lmm99 (= EMM99/tMMm) for polytropic en- 
velope structures; the light travel time or limb darken- 
ing corrections are not included. In order to compare 
our results with the semi-analytical solutions, we extract 
four following characteristics of shock breakout from the 
corrected and uncorrected models; Tc^peak: color tem- 
perature at i = 0, <imag: days until bolometric magni- 
tude declines by 1 magnitude after the bolometric peak, 
£'rad,imag: radiation energy emitted from t = — timag 
to t = iimag, and ipoak^ pcak bolometric luminosity. 
The properties of our models are summarized in Tabled 
Their dependencies on i?preSN and E are shown in Fig- 
ures [6^-[6)i and compared with the semi-analytic solu- 
tions. Since the semi-analytic solutions slightly depend 
on Moj that is different in the numerical models by a fac- 
tor of ~ 2. Possible ranges of the semi-analytic solutions 
are shown in Figures [S^-lBJi. 

Comparing the uncorrected models and the semi- 
analytic solutions, they are quantitatively different but 
the dependencies are roughly consistent. The models 
with larger i?proSN have lower Tc,pcak, longer iimag, and 
higher E'rad.imag- The models with higher E have higher 
J^cpcak, £^rad,imag, and Lpcak and shorter timag- The 
models with different Z are distributed along a sequence 
of the models with different i?proSN- This indicates that 
the variations with Z can be interpreted by the varia- 
tion with i?proSN and that the metallicity alters shock 
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1.5 0.5 1 1.5 0.5 1 1.5 

Days since bolometric peak (observer frame) [Days] 

Fig. 7. — Apparent g'-band LCs of the models at z = 0.2 (red), z = 0.5 (green), z = 1 (blue), z = 1.5 (magenta), z = 2 (cyan), z = 2.5 
(yellow), and z = 3 (black). No extinction and no IGM absorption are assumed. The panels are the same as in Figure[3] The horizontal 
line shows a 5cr detection limit in g' band for Subaru/Suprime-Cam 1 hour integration (gray, http://www.naoj.org/cgi-bin/spcam_tmp.cgi, 
assuming 0'.'7 seeing, I'.'S aperture, and 3 days from New Moon). 



breakout mainly through the variation of stellar struc- 
ture. Accordingly, the Z — 0.02 models could be applied 
even for shock breakout in stars with different Z if they 
have the same i?proSN-"'^^ 

The difference between Tc^peak and Tmm99 partly stems 
from the fact that the color temperature is different from 
the radiation temperature by definition when the opac- 
ity depends on frequencies and the SED is not blackbody. 
However, it is notable that the semi-analytic solution is 
roughly consistent with the numerical models if Tmm99 
is reduced by a factor of 1.5 (Figs.[S^ and[Bl3). The semi- 
analytic solutions give slightly higher values also for the 
other properties. They are in agreement with the numer- 
ical models if <mm99, ^'MM99, and Lmm99 are reduced by a 
factor of 2, 3, and 1.5, respectively (Figs.[nt-|B)i). The re- 
duction factors make it possible to approximately derive 

^® It is obviously better to perform a radiation hydrodynamics 
calculation for an evolutionary presupernova model with proper Z. 



progenitor properties from a direct comparison between 
observations and the semi-analytic solutions. Moreover, 
the qualitative consistency with the semi-analytic solu- 
tions supports the reliability of our numerical results. 

The light travel time and limb darkening corrections 
slightly reduce Tc^pcak and enhance i?rad,imag but does 
not change their dependencies. On the other hand, 
the corrections considerably change the dependencies of 
iimag on i?prcSN and E, and the dependence of Lpcak 
on E. This is because the corrections smear the LC 
peak and redistribute radiation energy emitted at bright 
epochs to a time range of Rph/c, where i?ph is a pho- 
tospheric radius. Consequently, the corrections lengthen 
iimag and diminish Lpeak more efficiently for a model 
with shorter timag and brighter ipcak- Shock breakout in 
the model with smaller -RprcSN and higher E has shorter 
iimag and brighter Lpeak (Figs.|ni;-[SJi and[S^-|nh) and thus 
more strongly corrected. 
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Fig. 8. — Apparent F322W-band LCs of the models at z = 0.5 (red), z = I (green), z = 
= 2.5 {cyan), and z = 3 (yellow). No extinction and no IGM absorption are assumed. 



1.5 (blue), z = 2 (magenta). 
The panels are the same as 



in Figure |3] The horizontal line shows a lOcr detection limit in F322W band for JWST/NIRCam 10* sec integration (gray, 
http: / /www. stsci.edu/jwst/instruments/nircam/sensitivity/index_html). 



3.2. Shock breakout at high redshift 

The theoretical multigroup LCs allow us to predict ap- 
parent LCs of SNe taking place at arbitrary distance, di- 
rection, and host galaxy. Let us consider the case where 
we observe an object at redshift z which emits light with 
a luminosity per unit frequency at frequency Pq in the 
rest frame and we detect the object with a flux per unit 
frequency /i/^to &t frequency Vohs in the observer frame. 
Here Vohs = ^'o/(l + z). The observed flux is obtained as 
follows: 



1 



obs 



(1 + z)L,„ (i/o) 



(1) 



where Dl is the luminosity distance. Here, we adopt 
for cosmological parameters a five -year result of Wilkin - 
son Microwave Anisotropy Probe (jKomatsu et al.ll2009l ): 
Ho = 70.5 km s"! Mpc-\ k = 0, Qx = 0.726, and 
r^M — 0.274. Apparent multicolor LCs are derived by 



convolving the diluted and redshifted multigroup LCs 
with the bandpasses of the satellites or telescopes. 

The LCs in g' band for z = 0.2, 0.5, 1, 1.5, 2, 2.5, 
and 3, LCs in F322W band for z = 0.5, 1, 1.5, 2, 2.5, 
and 3, and LCs in NUV and J bands for z = 0.05, 0.1, 
0.15, 0.2, 0.25, and 0.3 are shown in Figures [Tji-lThi, Ui- 
[5)ii,[5^-[5Jn, andfTUb-fTUtn. respectively. These figures also 
show the limiting magnitudes of telescopes/instruments 
with wide-field imaging capability: GALEX satellite 
in NUV band for 5a detec tion with 1500 sec integra- 
tion (rTiNuv.iim = 22.7 mag. iMorrissev et al.l[2005L I2007L 
Figs. [IK-On), Subaru/Suprime-Cam in g' band for 5(t 
detection with 1 hou r integration {nig' lim = 27.5 mag,^^ 
IMivazaki et allllM Figs. [7]a-l7ln), VLT/HAWK-I in J 



The limiting 
Subaru Imaging 



magnitude is calculated with 
Exposure Time Calculator 

(http://www.naoj.org/cgi-bin/spcam_tmp.cgi) assuming 0'.'7 
seeing, 1'.'5 aperture, and 3' days from New Moon. 
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Days since bolometric peak (observer frame) [Days] 

Fig. 9. — Apparent NUV-band LCs of the models at 2 = 0.05 (red), z = 0.1 (green), z = 0.15 (blue), z = 0.2 (magenta), z = 0.25 (cyan), 
and z = 0.3 (yellow). No extinction and no IGM absorption are assumed. The panels are the same as in Figure [3l The horizontal line 
shows a 5tT detection limit in NUV band for G ALEX satellite 1500s integration (ora^. Morris sev et al.|[2007l ). 



band for 5a detection with 1 hour inteKration (mjjim = 
24.8 mag. iKissler-Patig et al]l200 8l Figs. HOji-lMnj, and 
JWST/]<leaT Infrared Camera (NIRCam) in F322W band 
for 10a detection with 10''' sec integration (9.18 nJy,^^ 
Figs.lHK-lHlii). No extinction and no intergalactic medium 
(IGM) absorption are adopted here. 

These figures demonstrate that the shock breakout 
with higher E or larger i?proSN, i-e., larger Mms or higher 
Z, can be detected at higher redshift. The shock break- 
out can be detected in g' band even at z ^ 1 (13 and 
15 Mq models), z - 2 (18 and 20 M© models), and z ~ 3 
(25, 30, and 40 Mq). On the other hand, the observa- 
tions in NUV and near infrared (NIR) bands can detect 
shock breakout only at z < 0.5, although SNe at the later 
epoch, i.e., plateau stage, are detectable even at z > 4 by 
JWST (N. Tominaga, et al. in prep.). This is because the 
limiting magnitude in UV bands is much shallower than 

^^ |http: //www. stscl ■ edu/jwst/instruments/nlrcam/sensltlvlty/ 



that in optical bands and the SED of shock breakout 
is too blue for the NIR observations. The shock break- 
out has a blue color in optical (Figs. ITTk-fTllTL) . hence 
an observation in bluer optical bands is more suitable to 
detect shock breakout as long as the IGM absorption is 
irrelevant in the adopted bandpass. 

Shock breakout has large negative iC-corrections be- 
tween rest-frame x band and observer-frame x band , 
where x band is an arbitrary bandpass at A > 100 A. 
Figure fT2k shows i^FUV, ^nuv, Ku', Kg>, Kr', and Kj 
of the 20Mq, Z = 0.02, and S51 = 1 model at t = 0. 
Since the SED at t = is extremely blue, the negative 
i^-corrections are larger for higher redshift even in FUV 
band and the iiT-corrections in redder bands are smaller 
than those in bluer bands. Figure [T2b shows Kg/ of the 
Z = 0.02 models with different Mms at t = 0. The 
more massive models have slightly larger if-corrections 
by < 0.5 mag at z < 4.5. Figure fl2h shows the evolution 
in9t}!iKiijtiiQf the I5M0 and 3OM0 models. The if-correction 



10 




1.5 0.5 1 1.5 0.5 1 1.5 

Days since bolometric peak (observer frame) [Days] 

Fig. 10. — Apparent J-band LCs of the models. No extinction and no IGM absorption are assumed. The panels are the same as in 
Figure [3] and the color of lines are the sam e as in Figure (9] The horizontal line shows a 5cr detection limit in J band for VLT/HAWK-I 1 
hour integration {gray, Kisslcr- Patig et al.l l2008f). 



evolves more rapidly for higher redshift or smaller Mms- 
As a result, Kgt of the I^Mq model is smaller at t = 
but larger at iobs = 1-5 days than that of the 30Mq 
model, where tohs is a time from t = in the observer 
frame. 

Figures fT3k-IT3hi show distance modulus and apparent 
peak g'-band magnitudes mg'^poak of models as a function 
of redshift for different assumptio ns on the hos t galaxy 
extinction and IGM absorption (jMadaul I1995() . Here, 
we assume that the host galaxy has a color excess as 
E(B — V)host = or 0.1 mag and our Galaxy extinction 
law (IPeilll992D . These fi gures display that the dimming 
of apparent magnitude is considerably weak compared 
to the distance modulus because of the large negative 
ii'-correction. For illustration purpose, a 1-hour limit- 
ing magnitude for Subaru/Suprime-Cam is also shown 
('Tig'^iim = 27.5 mag). The maximum redshift for detect- 
ing shock breakout with 8m-class telescopes is mainly de- 
termined by the host galaxy extinction, while the IGM 
absorption becomes relevant at z > 2 in 5' band. When a 



limiting magnitude TOg'.iim ~ 30 mag is achieved by next- 
generation 30m-class telescopes, the maximum redshifts 
are z ^ 3.5 for E{B — l^)host = 0.1 mag and z ~ 4.2 for 
E{B — T^)host = which are limited by the IGM absorp- 
tion. At that era, shock breakout can be used for studies 
on not only distant SNe II-P, host galaxy extinction, or 
evolution of universe but also an IGM. 

3.3. Future prospects for shock breakout surveys 

The detection of shock breakout has recently been con- 
summated but the number is still small. This is because 
it is difhcult for past /ongoing SN/transient surveys with 
several days intervals to detect short-term soft X-ray/UV 
flashes. Hence, we propose a deep optical survey with 
short intervals, e.g., ~ hour, for shock breakout in an 
SN II-P in the distant universe. 

This is motivated by the following two brand-new 
prospects. (1) Due to the large negative X-correction, 
distant shock breakout up to 2: 3 is bright enough to 
be detected with current optical facilities 13. 2p . Such 
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Fig. 11. — Color {g' — r') evolution of models in the observer frame at z = 0.2 (redV z 
= 2 (cyan), z = 2.5 (black), and 2 = 3 (gray). The panels are the same as in Figure|3l 



0.5 (green), z = I (blue). 



1.5 (magenta). 



a survey is promising because of three reasons: the du- 
ration of a distant event is elongated, the star forma- 
tion rate (SFR) is high in the distant universe (e.g., 
iHopkins fc BeacorDll2006[ ). and available optical facilities 
are numerous compared to X-ray/UV satellites. (2) The 
multiepoch imaging observation in a night is essential to 
draw the LCs of shock breakout in both rising and de- 
clining phases because the time scale of shock breakout 
in an SN II-P at z < 3 are less than ~ 1 day in the 
observer frame 13.21) . 

In the following, we estimate the expected number and 
highest redshift of detection, discuss influences of uncer- 
tainties on host galaxy extinction and SFH, and propose 
realistic and promising survey strategies, ways to identify 
shock breakout, and ways to constrain the SN properties 
from observable quantities. In this section, we focus on 
the models with E^i = 1 and Z = 0.02. 

3.3.1. Expected number 



Host galaxy extinction heavily reduces brightness of 
shock breakout (§ 13. 2p . and thus it should be taken into 
account for a realistic number estimate. However, the 
host galaxy extinction of distant SNe II-P is unknown 
but instead it will be clarified by future shock breakout 
studies. Hence we expediently assume that the distri- 
bution of host galaxy extinction of distant SNe II-P is 
equivalent to that of nearby SNe II-P. The host galaxy 
extinction of a nearby SN II-P is estimated from Na I- 
D lines of the host galaxy, a spect roscopic observation of 
SNe, or a color of SN pla teau (e.g. JKrisciunas et al.ll2009l : 
lOlivares E. et al.l |2010[) . We employ the distribution 
of hos t galaxy extinction presented in lOlivares E. et all 
l|20Tnh . Although we take conservatively the highest ex- 
tinction for each SN among their estimates, we caution 
that it could be biased towards bright, i.e., less-reddened, 
SNe II-P. 

We estimate the incidence rate of shock breakout that 
can be brighter than a limiting magnitude mx,iim in 
bandpass x, named "observable SN rate" . The procedure 
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Redshift z Redshift z Days since bolomeric peak (observer frame) [Days] 

Fig. 12. — if-corrections Kx between rest-frame x band and observer-frame x band, (a) -ffpuv {I'^d,), i^Nuv (5''66")i (blue), Kgi 

(magenta), K^i {cyan), and i^j (black) of the 20Mq, Z = 0.02, and £51 = 1 model as a function of redshift. (b) Kgi of the Z = 0.02 and 
E^l = 1 models with Mms = 13Mq (red), 15Mq (green), ISM© (blue), 20Mq (magenta), 25Mq (cyan), 30Mq (black), and 40Mq (gray) 
as a function of redshift. (c) Kgi of the Z = 0.02 and £51 = 1 models with Mms = 15Mq (thin lines) and 30Mq (thick lines) at z = 0.2 
(red), 0.5 (green), 1 (blue), 1.5 (magenta), 2 (cyan), 2.5 (black), and 3 (gray) as a function of time in the observer frame. 
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Fig. 13. — Peak g' magnitude "ig'^poak s-^id distance modulus (black) as a function of redshift. The peak magnitudes are derived by 
assum ing no extinction and I GM absorptio n (red), host galaxy extinction (E(B — y)host = 0.1 mag and our Galaxy extinction law. IPeil 
119921) (blue), IGM absorption l lMadaulll995l ) (green), and host galaxy extinction (E(B — V)host = 0.1 mag and our Galaxy extinction law) 
and IGM absorption (magenta). The panels are the same as in Figure [3] and the horizontal line (gray) is the same as in Figure [T] 
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Fig. 14. — Observable SN rates per square degree per day 
nri;(ma;^liin) for (a) different bandpass: FUV {red), NUV {green), 
u' {blue), g' {magenta), r' {cyan), i' {yellow), z' {black), J {or- 
ange), and F322W {gray) and for (b) different host galaxy extinc- 
tion and SFH in g' band: estimates HB06-Gal {red), HB06-LMC 
{cyan), HB06-SMC {gray), MDP07-no {green), MDP07-Gal {blue), 
and MKll-Gal {magenta). The insets enlarge the regions with 
m^^jim < 28 mag. 

to estimate observable SN rate is as follows: we define 
a peak magnitude in bandpass x [TOpcak,x(-^MS) ^Vj -z)] 
for an SN with Mms which explodes in a galaxy 
with extinction Ay at a redshift z, and then a de- 
tection probability fx[mpcsk,x{MMS„ ^v, z),mx,i\m\ is set 
to 1 when mpoak,a;(MMS, ^v, 2) < rnx,\im or when 
mpcak,x(MMS,^v,2;) > rnx,\im- Accordingly, the observ- 
able SN rate per unit solid angle per unit time in the 
observer frame nx{'mx,\i-m) is obtained by integrating the 
detection probability with an IMF 0(Mms): a cosmic 
SFH ry(z), and the distribution of host galaxy extinction 
x{Ay), where j x{Av)dAY = 1, as 



fx [^pcak,a: ^) ; ^a;,lim] 



dV{z) r]{z) 
dQdz 1 -\- z 



x{Ay)dMMsdzdAy,{2) 



where V{z) is a comoving volume up to z and 57 is a 
solid angle. Here, we assume a modified Sa lpeter A IMF 
(lSal2eter"1955': Baldrv & Glazebrook"2003) and a cosmic 
SFH in Hopkins & Beacom (2006), and we adopt approx - 
imate formulae for the IGM absorption (jMadaulll995| ) 
and the host galaxy extinction (our Galaxy extinction 
law. lPdl[T992h .^9 

Since these formulae cannot be adopted at A < 800 A in 
the rest frame, we assume that the absorption and extinction at 
A < 800 A are the same as those at A = 800 A. 
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Fig. 15. — Reachable redshifts ^i,/ ("i^ nm) for (a) different band- 
pass with / = 0.9: the colors are the same as in Fig. 114b . for (b) 
different / in g' band: / = 0.999 {red), 0.99 {green), 0.9 {blue), 0.5 
{magenta) , and 0.1 {cyan), and for (c) different host galaxy extinc- 
tion and SFH in g' band; the colors are the same as in Fig. 114b . 
The insets enlarge the regions with 1;,,^ < 24 mag. 

The observable SN rate per square degree per day 
are shown as a function of mx.iira in Figures 114b . and 
[Mb . For deeper limiting magnitude, the observable SN 
rates are higher due to the larger observable volume and 
the higher SFR at high redshift. For example, the g'- 
band observations with mg/jim = 20, 22, 24, 26, 27, 28, 
and 30 mag detect 3.1 x 10~^ 6.2 x 10"^, 1.4 x 10"^ 
3.7 x 10-\ 1.7, 5.6, and 1.8 x 10 SNe degree-^ day-^, 
respectively. ngi{mgi \[^) logarithmically increases with 
Tig', Urn at TOg'.iim < 27.5 mag, but the increase slows 
down at m^'jim ^ 27.5 mag because the maximum red- 
shift of observable SNe reaches z ~ 2.5 (Figs. [T5k-fT5hi) 
where the IGM absorption becomes nonnegligible and 
the cosmic SFH hits a peak. 

Redshifts below which a given fraction / of observ- 
able SNe locates [zx,f(jnx,\i-m), named "reachable red- 
shift"] are shown in Figures [TSk-ITSh. The fraction of 
high-z events increases with mx,\ini- For example, the 
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Fig. 16. — Observable SN rate per day f^obs'^a:('"i:,lim) as a 
function of limiting magnitude m^, and survey area Cobs- (a) 
^obs'^a:("^2! lim) = 1 in different bandpass: the colors are the 
same as in Fig. I14h . and (b) different nobs"a;("^i!,lim) in g' band: 
^ohafixin^xfim) = 0.01 (red), 0.1 {green), 1 (blue), 10 (magenta), 
100 (cyan), 10^ (yellow), 10"* (black), and 10^ (orange). The lines 
being accessible by equal survey powers are shown (gray). 

g'-band observations with mg/jim = 20, 22, 24, 26, 27, 
28, and 30 mag result in [zg'.o.i(™g'.iim), V,o.5(™sMim): 
V 9(^3' lim)] = (0.013, 0.027, 0.042), (0.035, 0.071, 
O.l'l), (0.095, 0.19, 0.31), (0.26, 0.55, 1.0), (0.42, 0.96, 
1.9), (0.62, 1.4, 2.4), and (0.90, 1.9, 3.0), respectively. If 
?7ig', lim > 27 mag is realized, more than half of the ob- 
servable SNe take place at z > 0.9. This capability to 
access high-z universe is an intriguing and unique feature 
of shock breakout. 

If the same limiting magnitude is available, nx{mx,iim) 
and Zx,f(mx,iim) in bluer bands are higher than those in 
redder bands at shallow mx,\im, e.g., m^^jim < 24—25 mag 
for UV, due to the blue SED of shock breakout. However, 
the increases in nx{mx,iim) and Zx,f{'mx.\im) for deep 
Wx,iim are suppressed by the IGM absorption. Therefore, 
nxiiTix,!!!!!) and Zx,f{rnx,iiin) in redder bands overcome 
those in bluer bands for deep m^jjim. Thus, the most ef- 
fective bandpass for shock breakout detection depends 
on feasible mx,iiui- For example, the g'- and r'-band 
observations are currently the most effective in num- 
ber and in reachable redshift, respectively, because only 
8m-class optical telescopes or Hubble Space Telescope 
('Tlx. lim ^ 28 mag) are available at the moment. When 
30m-class optical/infrared telescopes {nix.nm ~ 30 mag) 
become operative, the most effective bands will be r' 
band in number and i' band in reachable redshift. 

The reachable redshift dramatically increases if mx,rnn 
is deeper than 26 — 30 mag (Figs. [T5k-IT5b) because 
the large negative iiT-correction makes the apparent 
peak magnitudes of shock breakout almost constant 



26 — 30 mag) for a wide redshift range (Figs. [T5k-ll"5hi) . 
Therefore, in order to detect distant shock breakout at 
z > 1, it is essential to attain m^^^um > 26—30 mag. Since 
the dramatical increases of reachable redshift in optical 
bands coincide with the limiting magnitudes of current 
optical facilities, the improvement in near future will en- 
hance the reachable redshift considerably. On the other 
hand, the increase in Za;,o.9(''TT'a;ji„i) reaches the ceiling at 
z ~ 4 for mx,iim > 30 mag even in NIR bands. This 
is because the cosmic SFH inESik ms fc BeacomI ([20061 ) 
has low SFR at such a high redshift. It is important to 
note that the detection of shock breakout at z > 4 is 
feasible if an SFR is high enough, unless the bandpass 
are below the rest-frame Lya wavelength. 

According to nx{mx,iim), observable SN rates per unit 
time i^obs'T-a;(TOx,iim) with a given mx.um and a survey 
area riobs are shown in Figures 116b . and 116b . An obser- 
vation in a bluer band is more efficient with the same 
'm,x,iiin if a wide survey area is available (> 1 degree^), 
while an observation in a redder band is slightly better 
if only a narrow survey area observation (< 0.1 degree^) 
with deep mx.nm (> 30 mag) is available. The fiat depen- 
dence of ilohsnx{mx,iim) on mx,iim at deep mx,iim, e.g., 
Wg'jim > 28 mag in g' band, stems from the suppression 
of nx{mx,iim) by the IGM absorption and the low SFR 
at high redshift. 

Figure fTBb shows lines giving equal ilobs'T-g' (™g',iim) 
with the g'-band observation. The gray lines represent 
equal survey powers without taking into account ob- 
serving overhead such as readout time. The number 
of observable SNe is larger for wider and shallower ob- 
servations with a given survey power. However, such 
a wide and shallow observation misses high-z events 
(Figs. [TSk-ITSb) and Hobs^ig' (wg'.iim) has an upper limit 
somewhere if the overhead is taken into account. There- 
fore, the practical survey parameters should be cus- 
tomized to purposes of observations and adopted tele- 
scopes/instruments, considering the number, reachable 
redshift, and overhead. 

3.3.2. Dependencies on host galaxy extinction and star 
formation history 

We set the ab o ve e stimate with the SFH in 
ms fc BeacomI (|2006[ 1 and our Galaxy extinction 
law (HB06-Gal) as a control estimate and investigate the 
dependence of rix{mx,iini) on uncertainties of host galaxy 
extinction and SFH fFig.lTIb). 

We attempt Large Magellanic Cloud and Small Mag- 
ellanic Clou d extinctio n laws (estimates HB06-LMC and 
HB06-SMC,|Pel[i99l) for host galaxies. The LMC and 
SMC extinction laws have larger absorption at A < 
2000 A and smaller at A ^ 2200 A in the rest frame than 
our Galaxy extinction law. As a result, Ugi {mgi um) of es- 
timates HB06-LMC and HB06-SMC are slightly smaller 
at TTig' Jim > 27 mag than that of estimate HB06-Gal but 
the estimates HB06-Gal, HB06-LMC, and HB06-SMC 
are similar at TOg'^um < 26 mag. 

The SFH in iHopkins fc BeacomI (|2006l ) is derived by 
scaling UV SFH so as to be consistent with infrared 
SFH and thus presumably includes both visible and 
dust-obscured star formation. Although we correct the 
dust extinction in host galaxies, it could underestimate 
the host galaxy extinction. Additionally, some studies 
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TABLE 3 
Simulated survey strategies. 



Strategy 


-^^night 


^obs 


^cxp 




— 

Number 


Reachable redshiff^ 


— 

Number with overhead 




[nights] 


[night- 1] 


[min] 


[mag] 


[degree--^] 




[degree-^] 


A 

A 


1 


3 


120 


28.4 


1.56 


2.40 


1.42 


B 


2 


3 


60 


28.0 


2.62 


2.21 


2.38 


C 


4 


3 


30 


27.6 


3.57 


1.99 


3.25 


D 


12 


3 


10 


27.1 


4.66 


1.46 


4.23 


E 


24 


3 


5 


26.7 


5.22 


1.07 


4.75 


F 


40 


3 


3 


26.4 


5.49 


0.84 


4.70 


G 


60 


3 


2 


26.2 


5.81 


0.71 


4.65 


H 


2 


4 


45 


27.9 


2.71 


2.25 


2.46 


I 


2 


5 


36 


27.8 


2.73 


2.25 


2.49 


J 


2 


6 


30 


27.6 


2.71 


2.25 


2.46 


K 


2 


12 


15 


27.3 


2.31 


2.12 


2.10 


L 


3 


1 


120 


28.4 


0.0651 


2.44 


0.0592 


M 


2 


3 


120 


28.4 


3.91 


2.35 


3.55 


N 


3 


3 


120 


28.4 


6.26 


2.34 


5.70 



^ The limiting magnitude is calculated with Subaru Imaging Exposure Time Calculator 
| |http : / /www ■ nao j . org/cgi-bin/sp cam_tmp . cgi[ ( assuming 0"7 seeing, aperture, and 3 days from New Moon. 

The number of observable SNc in total, assuming no overhead. 
^ The rcdshift below which 90% of observable SNe occur. 

^ The number of observable SNe reduced in proportion to the overhead, 10% for tcxp ^ 5 min and 30 sec for t^xp < 5 min. 



suggest that IHopkins fc BeacomI (j2006[) overcorrects the 
dust extinction and overesti mates the SFR by a factor 
of 2 at z 2 - 3 r e.g.. iNagashima fc Yoshiil 120041: 
Nagashima et alJ 120051: iBaugh et alJ 120051: iLacev et alJ 



20101) ■ Hence, we tes t with UV dust-unobscured SFH 



( Mannucci et al.ll2007[) . Although their estimate is hni- 
ited at z < 2, the SFR at z > 2 is assumed to be the same 
as the SFR at z = 2. Since the dust attenuation of shock 
breakout brightness in a host galaxy would be somehow 
involved in the UV dust-unobscured SFH, we present an 
estimate with no additional host galaxy extinction (esti- 
mate MDP07-no) and also test an estimate with extinc- 
tion correction with our Galaxy extinction law as a lower 
hmit (estimate ]VIDP07-Gal). 

Furthermore, we test an SFH and host galaxy extinc- 
tion derived from a cold dark matter-based semi-analytic 
mod el (M. A.R. K obayash i, et al. 2011 in prep., see 
also INagashima k Yoshii 2 00l iKobavashi et al.l 120071 
12010)'^'^ and our Galaxy extinction law (estimate IVIKll- 
Gal). The estimate MKll-Gal provides a self-consistent 
observable SN rate with respect to the host galaxy ex- 
tinction and SFR, although the metallicity evolution is 
ignored in this paper. The estimates HB06-Gal, 1VIDP07- 
no, and MKll-Gal are consistent but they are higher 
by a factor of 2.5 than the estimate ]V[DP07-Gal at 
"T-g'.iim < 28 mag. On the other hand, the differences be- 
tween the estimate HB06-Gal and the estimates MKll- 
Gal and MDP07-no are as large as a factor of 2 — 2.5 
at TOg'jim ~ 30 mag. Therefore the number count of 
shock breakout with 30m-class telescope can constrain 
how high the SFR is at high redshift, being independent 
of galaxy studies. 

We also investigate the dependencies of Zg>^ f{m n',iin i) 
on host galaxy extinction and SFH (Fig. Il5b). 
^g',o.9("^g',iim) of the estimates HB06-Gal and MDP07- 
Gal are consistent but they are slightly lower than the 

The model well explains observations of nearby and distant 
galaxy evolutions and dust-unobscured luminosity density. Since 
the intrinsic SFR and dust extinction of each galaxy are provided, 
the SFH estimate is free from uncertainties of extinction correction 
and conversion from the galaxy luminosity to the SFR. 



estimate MDP07-no at mg/.iim > 26 mag and the es- 
timate MKll-Gal at m^'jim < 25 mag. And the esti- 
mate HB06-LMC has slightly higher Zg',o.9('7ig',iim) at 
'7^g',lim < 25 mag than the estimate HB06-Gal, while 



,o.9("T-g',iim) of the estimate HB06-LMC at m 



> 



27 mag and the estimate HB06-SMC are consistent with 
that of the estimate HB06-Gal. These are because 
Zxji'Tix.iini) depends mainly on the host galaxy extinc- 
tion at shallow mxAim and the IGM absorption at deep 
'Tiajjim- The reason why Z£!;,o.9('^a;,iim) of the estimate 
MDP07-Gal is higher than that of the estimate HB06-Gal 
at mx^iim > 31 mag is because we assume the constant 
SFR at high redshift in the estimate MD P07-Gal wh ich 
is higher than the SFR in , Hopkins fc BeacomI (|2006[) at 
z > 5. 



3.3.3. Requirements on survey strategies 

Although the above estimates refer only to the peak ap- 
parent magnitude, the cadence of observations is also an 
important ingredient to identify transients. Indeed, it is 
difficult to identify and interpret an event only with one- 
epoch brightening. Additionally, the above estimates do 
not take into account an elongation of the duration at 
high redshift. Hence, we estimate the number of detec- 
tion of shock breakout Ni^^ {iTix.iiim^obs) for an observa- 
tion sampling i with bandpass x, mx^iim, and a field of 
view r^obs, which can be obtained by an integration of a 
"control time" Ti^xiMus, ^v, z, nix, lim)^^ as follows: 



A^z.a: (^a;,lim; ^obs) 



r,,^(MMs,^v,2 



dV(z) ri(z) 

0(MMs)»obs \'\' x{Av)dMMsdzdAy,{3) 
ailaz 1 + z 

which is similar to Eq. ([2]) but integrates 
with Fi ,j.(AfMSi ^Vj ^7 ™a;,iim) instead of 

Fi^j; (MmSi ^Vi -^i "lai.lim) is the number of detection of shock 
breakout with MmSi ^Vi ^^iid ^ per unit SN rate in the observer 
frame and thus having a dimension of number/(number/time), i.e., 
time. 
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fx[mpcak,x{MMS,Av,z),ma:,iim]- We define the de- 
tection with two criteria: (1) the event is detected at 
> A^dctcct samplings with 3cr and (2) at least one of sam- 
plings is taken from tobs = —0.2 days to ^obs = 0.4 days. 
ri,x(-^MSi ^V: -z, m^, is determined by simulating 
the detection of shock breakout. The expected number 
of detection and reachable redshift are slightly smaller 
and lower than the previous estimates in which all 
events with TOpeak,a: < iTT-xMrn are counted because we 
require the detectable SN to be brighter than m^jjim 
over > iVdctcct samplings. 

We attempt strategies A-N which divide a 6-, 12-, or 
18-hour (/'-band observation with respective ways, as - 
suming the cosmic SFH in IHopkins fc BeacomI (gOOl) 
and our Galaxy extinction law ()Pei| 119921) . Here, we 
treat the number of nights A'night, observation sam- 
pling per night nobs, and exposure time for each sam- 
pling ^exp as parameters, assume A'detect = 3, and adopt 
Suba ru /Suprime- Cam to e stimate nix^iim for 3cr detec- 
tion (iMivazaki et al.l [20020 . Each observation sampling 
is uniformly distributed to the whole night (10 hr for each 
night) and an overhead is neglected. The detail of strate- 
gies and the results, the expected number per square de- 
gree and the redshift below which 90% of events occur, 
are summarized in Table [31 

While large A'night enhances the number of detection, 
large riobs and long iexp enhance the reachable redshift. 
For the same A'night, a strategy with nobs = 5, in which 
a field is observed every 2 hours, is the most efficient in 
number and reachable redshift. The remarkably small 
number of detection of strategy L (nobs = 1) demon- 
strates that the multiple photometric observation in a 
night is essential to detect shock breakout because shock 
breakout brightens and declines within ~ 1 day. Com- 
paring the strategies A, M, and N, an additional one 
night observation increase the number of detection by 
2.3 degree"^, which is larger than the number of detec- 
tion of the one night observation (1.6 degree"^, strategy 
A). This is because the observation over > 2 continuous 
nights can detect an event, taking place at the end of the 
first night, at the subsequent nights. Thus, the strategy 
with iVnight > 2 is more favorable than the strategy with 

^night — 1- „^,„ 

As expected in § 13.3.11 a wider and shallower survey 
detects larger number of shock breakout at lower red- 
shift. However, the fraction of overhead should be high 
for extremely short icxp in reality. Thus, we estimate 
the number of observable SNe taking into account an 
overhead. Here, we expediently assume an overhead as 
10% of total observation for t^xp > 5 min or 30 sec for 
texp < 5 min, and reduce the number of observable SNe 
in proportion to the fraction of overhead. As a result, 
the number of observable SNe peaks at tcxp — 5 min for 
this specific overhead. 

3.3.4. Identification of shock breakout 

Many transients other than shock breakout, e.g., vari- 
able stars, SNe, quasars (QSOs), and GRBs, will be 
found as variable objects in a photometric observation. 
Shock breakout can be reliably discriminated from other 
kinds of variable objects holistically referring to observ- 
able quantities: time scale, LC shape, color, and position, 
and other observations. 
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Fig. 17. — Color-color {u' —g' vs. g' — r') diagrams for the models 
with Mms = 15M0, Z = 0.02, E51 = 1 (red), Mms = 25Mq, 
Z = 0.02, E51 = 1 {green), and Mms = 25M0, Z = 0.02, E51 = 10 
(blue), at (a) z = 0.5, (b) z = 1, (c) z = 1.5, (d) z = 2, (e) 
z = 2.5, and (f) z = 3, without the IGM absorption {thin line) 
and with the IGM absorption (thick line). The points represent 
the colors at t = {circle) and at t = 2 days in the observer frame 
{square). The contours show the distributions of stars {cyan, SDSS, 
Abazaiian et al. 2009) and QSOs {magenta, Schneider et al. 200^, 
2Q03, 2005, 2007,, ^010). The arrows show extinction vectors of 
host galaxies with Ay = 0.25 mag, assuming that the extinction 
curve is th e same as our Galaxy {red), LMC {green), or SMC {blue) 
UPeai 19921 1. 



Time scale and LC shape: Shock breakout has a 
nonrecurrent brightening, a time scale of several sec — 
several days, and featureless LC. On the other hand, 
(1) SNe powered by radioactive decays, the plateau of 
SNe II-P or Type Iln SNe, or linear decay of Type II-L 
SNe have time scales of several ten days — several hun- 
dred days, (2) GRB optical flash has a time scale of sev- 
eral t en sec and fre quently has a jagged multiple peaks 
(e.g.. lWozniak et al., 2009). and (3) the flare of low -mass 
star, e.g., a M dwarf star, has a time scale of several 
minutes — several h ours and recurrent brightening (e.g., 
iHawlev et al.l [20031) . Therefore, the time scale and LC 
shape of a transient can be employed for excluding these 
objects. 
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g' - r' [mag] 

Fig. 18. — Color-color ( g' — r' vs. r' — i') diagrams for the models at various redsliifts. The panels, contours, points, lines, and arrows 
are the same as in Figure [TTl but A-y = 0.5 mag is assumed for the extinction vectors of host galaxies. 




-0.5 



0.5 



-0.5 



0.5 



-0.5 



0.5 



FUV - NUV [mag] 



Fig. 19.— Color-color (FUV-NUV vs. g' - r') diagrams of the models at (a) z = 0.05, (b) z = 0.1, (c) z = 0.15, ( d) z = 0.2, (e) 
z = 0.25, and (f) z = 0.3. The points without li nes represent colors of sta rs (cyan, Bruzual-Persson-Gunn-Stryker atlas. IStrecker et al.l 
11979 : [Ciunn &: Strvkerl ll983') and QSOs ( magenta. Wanden Berk et al.|[20011) . The colors of QSOs are represented by open triangles w hen 
the rest-frame Ly a wavelength is redshifted into FUV band. The color of lines, points on lines, and arrows are the same as in Figure [17] 
but Ay = 0.5 mag is assumed for the extinction vectors of host galaxies. 

Position, archival image, and other observa- 
tions: If a survey is performed at fields with plenty past 
observations, checking past variabihties at the position 
can effectively rule out possibilities of long time-scale 
variables such as QSO or a variable star. If the event oc- 
curs in the outskirts of the host galaxy or is not detected 
in X-ray, the variable object is likely to be shock break- 
out. And the deep UV, optical, and infrared imaging 

Some low-luminosity AGN at high rcdshift can not be detected 
even in deep X-ray data (,Saraiedini ct al 200&; .Cohen et al 200^ : 



data are also useful for excluding the possibility of stars. 
Furthermore, if the field is include d in the field of view 
of 7-ray telescope (e.g., 5'wzft/BAT. lGehrels et aIl l2004V 
an alert of GRB can be used for ruling out a possibil- 
ity of GRB. Even if the GRB prompt emission cannot 
be observed (i.e., the GRB orphan afterglows), a radio 
follow-up observ ation can constrain the prese nce of rela- 
tivistic jets I'e.g- ISoderberg et al.ll2008ll20lB . 

IMorokuma et al.ll2008t) . 
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Fig. 20. — Color-col or ( J — H vs. H — K) diagrams of the models at various redshifts. The panels, color of lines, points, and arrows are 
the same as in Figure [191 but Ay = 3 mag is assumed for the extinction vectors of host galaxies. 



Color: The conclusive identification of sliock break- 
out is given by tfie color of a variable. Figures [TTk - 
[nt {u' - g' vs. g' - r'), {g' - r' vs. r' - 

i'), [Hi^ (FUV-NUV vs. g' - r'), and [20]a-l20t 
{J — H vs. H — K) show redshift-dependent color- 
color diagrams for the less- massive ISM© and Ec,i — 1 
model, massive 25M0 and i^si = 1 model, and ener- 
getic 25M0 and E^i = 10 model with t aking into ac- 
count the IGM absorption (jMadaul '1995') and compare 
them with colo r-color distrib utions of stars^^ a nd QSOs 
(jSchneider et a l. 2002, 2003!, 12001 l20?)7l IMoh for Fig- 
ures [TTk-fTTF and flSk-ITSt' and with colors derived from 
typical sp ectra (stars: Bruzual-Persson-Gunn-Stry ker at- 
las, Strcckcr ct al.' 197^; 'Gun n fc Str^^ |lp3p4 
QSOs: Vandcn Berk ct al. 2001) for Figures [ma-Ulf and 
I^Ok-Ent. Here, these figures show extinction-uncorrected 
colors of stars and extinction-corrected colors of QSOs 
at z ± Az [where Az — 0.1(1 -I- z), being an accuracy 
of a photometric redshift] . Ordinary SNe are not shown 
because any type of SNe are tog red to appear in these 
figures (e.g., INugent et all 120021) . The colors of shock 

breakout are FUV-NUV 0.5, u'-g' < -0.3, g' -r' < 

»■' -0.3, J-iJ < -0.6, andiZ-iiT < -0.8at 



t = and FUV 



NUV- 0.1, u'-g' < 0.1, g'-r' < -0.1, 
I' < -0.1, J ~ H < -0.3, and H - K < -0.4 
at tohs — 2 days (see also Figs. ITTb-fTThi). except for 
u' — at z > 2.5 or — r' at z > 3. These exceptions 
stem from the fact that the light in observed u' and g' 
bands are heavily absorbed by the IGM at z > 2.5 and 
z > 3, respectively. The color of shock breakout is much 
bluer than the majorities of stars and QSOs but we note 

Photometrically-identified stars with photometric errors < 
0.01 mag are extracted from the Seventh Data Release (DR7; 
|S?azaiian et al. 2009) of the Sloan Digital Sky Survey (SDSS; 
lYork et al.ll2000l ). The number of stars is over 4 X 10^. We note 
that the sample could be biased toward bright and thus blue stars 
because the stars with small errors are selected. 

•^"^ [http: //«w. stscl ■ edu/hst/observatory/cdbs/astronomlcal. 



that the color at t = is similar to 0/B stars and that 
the NIR colors of less-massive and energetic models at 
iobs = 2 days and massive model at tohs = 10 days are 
similar to QSOs. 

GRB orphan afterglow and M dwarf flare could have 
a similar time scale, featureless LC, and nonrecurrent 
brightening in a limited-time photometric observation. 
However, the blue color is a precise identifier of shock 
breakout. According to a standard model for GRB af- 
terglow (e.g., fSari ct al. 1998), GRB afterglow are red 
at a frequency range above frequencies corresponding to 
minimum or cooling Lorentz factor. Adopting reas onable 
parameters for GRBs (e.g.. iPanaitescu fc Kumail [2002). 
the minimum or cooling frequency at > 0.1 days after 
the prompt burst is v < 10^'^ Hz and thus the red color 
is realized at A < 3 x lO'' A. Although the minimum 
and cooling frequencies are higher for earlier epochs, the 
GRB orphan afterglow will typic ally peak at > 0.1 days 
(e.g.. iTotani fc Panaitescul l2'002( ) . Therefore, the GRB 
orphan afterglow is likely to have a red color in optical 
bands. Also the color of low-m ass star flare is typically 
red fe.g.. lKowalski et al.|[2009l) . 

3.3.5. Constraints on supernova properties from shock 
breakout 

The SN properties can be constrained from the ob- 
servations of shock breakout. Here, we focus on three 
observational quantities: peak magnitude, decline rate, 
and color evolution. Their dependencies on the model 
properties, Mmsi -RprcSN, and E, ai z = 0.2, z = 1, and 
z = 2 are summarized in Tables |4l[5l|6] and shown in Fig- 
ures[2Tk-l2Tb.[22k-l22l. and[23k-l23l. respectively. Here, we 
assume that the redshift of shock breakout is determined 
photometrically or spectroscopically. 

Figures l^-I^Tb.B^-l^. andl^-l^ show the appar- 
ent peak ^'-band magnitude mg/^pcak- Wg'^peak vary over 
_ ^1.5 — 2 mag depending on i?prcSN and thus Mms; al- 
f^t^Rgffr^S.n'tpgg^], are similar for the models with £'51 = 1 
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Fig. 21. — Dependencies of observational quantities on the model properties, Mms 
represents Z = 0.02 models (red), Mms = 20Mq models with Z 
E51 = 1, 4, 10, and 20 (green), (g-o) The colors of the models at t 
and 2 days {open circles). 



^preSNi El'iid E, at z = 0.2. The color of symbols 
0.001, 0.004, and 0.05 (blue), and Mms = 25Mq models with 
ohB^ = (filled circles), 0.5 days (open square), 1 day (open triangles) , 
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(Fig. |6^). On the other hand, although ipcak varies by 
an order of magnitude depending on E (Fig. \Ejp), the 
'7^3', peak range of the models with different E is only 
~ 0.7 mag. The different behavior between bolometric 
and monochromatic luminosities stems from the differ- 
ent photosphcric temperatures leading to different Tc. 
The models with larger i?proSN have lower temperature 
and thus redder SEDs and slightly brighter at A > 300A 
(Fig. 2^), while the models with higher E have higher 
temperature and thus bluer SEDs but the luminosities 
are almost similar at A > 800A (Fig. |3)d). We note that 
the g'-band luminosities peak at iobs ^ O.OOf — 0.05 days 
because of the temperature evolution (see also Figs. [7^- 

mn). 

Figures [HJl-EIf , [21i-112f , and [^Sl-iSF show days after 
the g'-band peak until gi'-band magnitude declines by 
0.5 magnitude tg',o.5mag- The dependencies of the de- 
cline rate in g' band on i?proSN (-^Ms) and E are almost 
similar to those in the bolometric LC. This is because the 
decline rate depends on energetics of SNe but not on the 



temperature when the SED peaks at shorter A than g' 
band. Since the dependencies of TOg',poak and tg',o.5mag 
on i?proSN (Mms) and E are different, in principal they 
could determine i?preSN (Mms) and E independently. 

However, apparent brightness is dimmed by extinction 
and thus having large uncertainties. Hence, we introduce 
the color evolution to resolve the uncertainties. Although 
the absolute color is also strongly reddened by extinction, 
the color evolution does not suffer from extinction unless 
extinction changes with time. Figures [2T^-l2Tb. [22^-1221. 
and[23|-l23l show the color evolutions, FUV-NUV, u'-g', 
and g' — r' for z = 0.2 and g' — r' for z = 1 and z = 2, 
respectively (see also Figs.fTTk-ITTln). These figures show 
the bluest color, that is realized at fj^"^ = correspond- 
ing to iobs ~ 0.03-0.2 days, and the colors at i™'""^ = 0.5, 
1, and 2 days, where color is FUV-NUV, u' — g' , or g' — r' . 
The colors of some models evolve toward red at first and 
then get back to blue again, which is shown as a loop 
structure in Figures [T3i-[nF, [HJi-lIHF, [HJi-lIlF, and W^- 

[M- 

The model with smaller i?proSN lias a bluer color at 
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Fig. 23.— Same as Figs. [21] but for z = 2. 



the peak because of the higher temperature and its color 
evolution is more rapid. The color evolution within 
2 days are A{g' — r') > 0.3 mag for the model with 

^prcSN < IO^Rq and A(g'-r') < 0.25 mag for the model 
with i?prcSN > lO'^i?0. For example, the ISMq model 
has a redder color than the 30Mq model at t^^°^ = 1 day 

for z = 0.2 and at t^^^^ = 2 days for z = 1 and z = 2. 
On the other hand, the bluest colors of models with dif- 
ferent E are similar but the models with higher E have 
more rapid and smaller variations. The color evolutions 
of the models with different Z are similar to those of 
the Z = 0.02 models with the same -RprcSN and thus it 
could cause an uncertainty for constraining Mms when 
the metallicity of host galaxy is unknown. In conclu- 
sion, the color evolutions can classify shock breakout, 
at least, to the following three groups: explosions with 
i?proSN < IO^Rq and i?p,.cSN > lO^i?© (Mms < 20Mq 
and Mms > 25Mo for the Z = 0.02 models) and an 
energetic explosion with E^i > 10. 



4. CONCLUSIONS & DISCUSSION 

Shock breakout is the brightest event in the SN with a 
Shockwave and could consummate the detection of CC- 
SNe in the high-z universe. We present multicolor LCs 
of shock breakout in SNe II-P with various Mms> Z, 
and E based on realistic stellar models. Using our theo- 
retical models, we investigate the dependencies of shock 
breakout properties on the progenitors and explosion en- 
ergies and present thorough prospects for future surveys 
of shock breakout. It is essential for identifying and in- 
terpreting shock breakout to observe a field more than 
once in a night in multiple blue optical bands, preferably 
over > 2 continuous nights. And, adopting standard cos- 
mic SFH, IMF, extinction distribution of host galaxies, 
the ^'-band observable SN rate for mg/^um = 27.5 mag is 
3.3 SNe degree"^ day~^ and the half of them locates at 
z > 1.2. 

We calculate 13 SN models with Mms = 13 - 4OM0, 
Z = 0.001 - 0.05, and E51 = 1 - 20 (Tab. [2|). The 
model with larger i?picSN, thus typically larger A^mS: has 
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lower Tcpoak, longer iimag, and higher £;rad,imag, while 
the model with higher E has higher Tc,poak, ^^rad.imag, 
and Lpcak- The metallicity affects shock breakout mainly 
through altering the stellar structure. The variations of 
Tcpcak, iimag, ^^rad, imag , and ipoak amoug the adopted 
models are - (2 - 5) x 10^ K, - 0.01 - 0.6 days, ~ (0.7 - 
5) X 10''* erg, and ~ (0.5 - 7) x lO'''^ erg s~^, respectively. 
The dependencies of numerical results are similar to 
those of the semi-analytic solutions (Matzner & McKee 
|1999() . The semi-analytic solutions are nearly consistent 
with the numerical results if TMMgg, ^'MM99, ^mm99, and 
iMM99 are reduced by factors of 1.5, 3, 2, and 1.5, re- 
spectively. 

According to the models, we predict the observational 
quantities of high-z shock breakout. Since shock break- 
out has a blue SED peaked at ~ 100 A, brightness 
of shock breakout at a fixed observed bandpass is less 
dimmed compared to the geometrical dilution, i.e., shock 
breakout has a large negative _fir-correction. This makes 
it possible to detect high-z shock breakout, e.g., up to 
z ~ 3.5 in g' band with 8m-class telescopes if there is no 
extinction in the host galaxy. Although shock breakout 
strongly suffers from the extinction and IGM absorption, 
it can be detected up to z ~ 2 if the host galaxy has the 
color excess E{B — F)host — 0.1 mag and our Galaxy 
extinction law. 

Convolving the Salpeter's IMF, cosmic SFH, host 
galaxy extinction, and IGM absorption, we estimate the 
observable SN rate as a function of bandpass and lim- 
iting magnitude. As a result, considering the opera- 
tive telescopes/instruments, the g' band observation is 
currently the most effect ive in the number of detection. 
Adopting cosmic SFH bv lHopkins fc BeacomI (2006) and 
our Galaxy extinction law, the observable SN rate is 
3.3 SNe degree"^ day~^ for mg'jim = 27.5 mag. Even 
taking into account uncertainties on host galaxy ex- 
tinction and SFH, the observable SN rate is > 0.93 
SNe degree"^ day~^. 

We also present the redshift distribution of observable 
SNe. For currently-available mg'^nm = 27.5 mag, 50% 
of observable SNe take place at z > 1.2. Furthermore, 
for mg'^iim = 30 mag, ~ 10% of observable SNe locate 
at z > 3. Since the reachable redshift increases dra- 
matically if TOa:jini > 26 — 30 mag is feasible, the next- 
generation telescopes/instruments will considerably en- 
hance the reachable redshift. The reachable redshift is 
almost independent of the uncertainties involved in host 
galaxy extinction and SFH. Therefore, the shock break- 
out is the most appropriate phenomenon to aim at the 
detection of high-z CCSNe. The first detection of normal 
CCSNe at z > 1 can certainly be achieved by the obser- 
vation of shock breakout in SNe II-P. The direct obser- 
vation of normal CCSNe at z > 1 will shed light on their 
nature and cosmic evolution histories, most of which are 
currently derived from galaxy studies that might be bi- 
ased by brightness of galaxies. 

Future/ongoing wide and / or deep surveys, e.g., Palo- 
mar T ransient Factory fPTF. lLaw et al.ll2009l: iRau et al.l 
2009f). Lic k Observatory Supernova Search (LOSS, 
Leaman et al. 201 Qi), Catalina R eal-Time Transient Sur- 
vey rCRTS. iDrake et al.l[2009l ). Kiso/Kiso Wide Field 



Camera (K WFC),^^ Skymapper,^ ^ Dark Energy Sur- 
vey (DES, iBernstein et al.l 120090 . Panoramic Survey 
Telescope a nd Rapid Response System (Pan-STARRS, 
I Kaiser et~a l. 2002; Kaiser 2004), Subaru/Hyper Suprime- 
Cam (HSC, Mivazaki ct al. 200 (j), and Lar ge Synoptic 
Survey Telescope (LSST, Ivezi c et al.l 120081 ) . will find a 
large number of shock breakout. We simulate realis- 
tic survey strategies and show that a wider and shal- 
lower survey leads to a higher observable SN rate with 
a given survey power but misses the higher-z events. 
And the observable SN rate for short integration is sup- 
pressed by the overhead. Although the survey parame- 
ters should be customized to observation purposes and 
telescope/instrument, we conclude that the most essen- 
tial observation is the multicolor photometry with short 
intervals less than 1 day and that the observation over 
> 2 continuous nights is favorable. 

We also establish the ways to identify shock breakout 
and to constrain the SN properties from the observa- 
tions of shock breakout. The LC, color, position, and 
past or wide frequency-coverage observations excellently 
distinguish shock breakout from the variable stars, SNe, 
GRBs, and QSOs. In particular, the blue color of shock 
breakout is the most important information to identify 
shock breakout. Shock breakout in an SN II-P with 
larger i?preSN evolves more slowly and has a more lu- 
minous peak in optical bandpass. The two observational 
quantities, time scale and color evolution, can reasonably 
determine the SN properties, -RproSN (-^^Ms) and E, be- 
ing independent of the host galaxy extinction. When nu- 
merous shock breakout is detected, the IMF in the high-z 
universe will be constrained by shock breakout. Further- 
more, if the SN properties can be determined only from 
the time scale and color variation with time, the other 
two observational quantities, peak magnitude and ab- 
solute color, can determine the host galaxy extinction. 
Combining the observations of the host galaxies, the re- 
lation between the host galaxy and its stellar contents 
might be constrained. 

Shock breakout enables an untargeted CCSN survey at 
unprecedentedly high redshift, which can provide large 
uniform CCSN samples with the comparable redshift 
range with GRBs. This would allow us to clarify the 
relation between GRBs and star formation and thus the 
GRB progenito rs. GRBs are hosted in blue, faint, and/or 
low-Z g alaxies ("Le Floc'h et al.ll2003tlkewlev et ai]l2007l : 
iLevesqu e ct al. 20j3) and thus have been suggested 
to require low-Z prog enitors. On the other hand, 
iKocevski fc WestI (|2010l ) and iMannucci eFall (|2010[ ) re- 
cently suggest that the GRB hosts follow a tight corre- 
lation among stellar mass, metallicity, and SFR of field 
galaxies and that the characteristics of GRB hosts can 
be explained only by the large SFR in \ow-Z galaxies. If 
the latter is correct, the host galaxies of CCSN should 
share the same properties as the GRB hosts because the 
lifetimes of GRB and CCSN progenitors are similar. 

http : //www. ioa. s .u-tokyo . ac . jp/kisohp/ top_e.html| 

http : //www .mso . anu. edu. au/skymapper/ 

For example, the lifetimes of stars with Mms ^ 4OM0, 
a putative GRB progenitor, and Mms ^ IOMq, a putative 
CCSN progenitor, are less than 5 Myr and 20 Myr, respectively 
(Schaller ct al. 1992). In order to make a difference between GRB 
and CCSN hosts, the star formation in the host galaxies should 
cease 5 — 20 Myr ago or start < 5 Myr ago. Although it is still un- 
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However, it is suggested that the GRB hosts are fainter 
and more irregular a t z < 1.2 and lower - ^ at z < 0.3 
than the CCSN hosts ()Fruchter et al.ll2006l:lModiaz et al.l 
l2008t). although there are rooms for improvement, e.g., 
the statistics and uniformity. Hence, the large CCSN 
sample of the same quality as GRBs, especially on sur- 
vey method and redshift range, could give unique infor- 
mation on the environment of GRBs. The comparison 
between CCSN and GRB hosts at z > 1 can provide 
an essential clue to unpuzzlc the issue. It will probe 
whether an SFH estimate with GRBs, which is extendible 
to z > 8, are biased or not. 

Aspherical shock breakout in a cocoon or a jet 
is sugg ested for Type Ic SN 2006ai/GRB 060218 
(e.g., Soderberg et al. 2006'; 'Camoana e t al.l 120061 : 
iWaxman et al. 2007; Ghisclhni et al. 2007) and Type lb 
SN 2008D/XR F 080109 (e.g., iSoderberg eFaLr[2Q08l: 
iMazzali et al.l [2008' ) . The asphericity stems from the 
compactness of the progenitor which leads to relativistic 
outflow at the tip of Shockwave. Therefore, in order to 
precisely deal with aspherical shock breakout, a multidi- 
mensional relativistic radiation hydrodynamics calcula- 
tion is required but has not been attained so far. On the 
other hand, polarization observations demonstrate that 
SNe H-P we focused have a spherical structure at the 
plateau phase even with an aspherical inne r core (e.g., 
[Leonard et all [20061 IWang fc Wheeleril200 8^ This is be- 
cause the progenitor of an SN H-P has a thick H envelope 
diluting the asphericity. Therefore, assuming the spher- 
ical symmetry is reasonable for SNe H-P and the results 
are applica ble for all SNe H-P. 

Recently iSmartt et al.l ()2009[ ) suggests that there is a 
maximum mass for the progenitor that can be SNe H- 
P [Mms < 16.5 ± 1.5Mq).'^^ This suggestion cautions 
the existence of a massive SNe H-P with Mms > 2OM0. 
However, the constraints are limited only to SNe H- 
P having occurred in nearby universe, in which indi- 
vidual stars can be resolved. It is not investigated 
whether such a low maximum mass for SNe H-P ex- 
ists in high-z universe or at low-Z environment. Also, 
the reason why the stars with Mms > 2OM0 cannot 
be SNe H-P is under debate. Some studies suggested 
that such a m assive star forms a black hole directly (e.g., 
iSmarttI I2OO9O or explodes as other kinds of SNe due to 
a strong mass loss and/or rotation (e.g., ISmartt et al.l 

der d ebate how short duration o f star formation activi ty is allowed 
(e.g., IMas-Hesse Hz Kunth|[T999l : IMcOuinn e t al."2010'). such short 
time scales are comparable with the age sp read s of OB associations 
and the lifetimes of giant molecular clouds IIMcKee fe Ostrikerl20071 
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Observable quantities at z = 0.2. 
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TABLE 5 
Observable quantities at z = 1. 
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TABLE 6 
Observable quantities at z = 2. 
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